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ABSTRACT:  In  situ  characterization  of  minute  amounts  of  fluids  that  rapidly  change 
their  rheological  properties  is  a  challenge.  In  this  paper,  the  rheological  properties  of 
fluids  were  evaluated  by  examining  the  behavior  of  magnetic  nanorods  in  a  rotating 
magnetic  field.  We  proposed  a  theory  describing  the  rotation  of  a  magnetic  nanorod  in 
a  fluid  when  its  viscosity  increases  with  time  exponentially  fast.  To  confirm  the  theory, 
we  studied  the  time-dependent  rheology  of  microdroplets  of  2-hydroxyethyl- 
methacrylate  (HEMA)/diethylene  glycol  dimethacylate  (DEGDMA) -based  hydrogel 
during  photopolymerization  synthesis.  We  demonstrated  that  magnetic  rotational  spectroscopy  provides  rich  physicochemical 
information  about  the  gelation  process.  The  method  allows  one  to  completely  specify  the  time-dependent  viscosity  by  directly 
measuring  characteristic  viscosity  and  characteristic  time.  Remarkably,  one  can  analyze  not  only  the  polymer  solution,  but  also 
the  suspension  enriched  with  the  gel  domains  being  formed.  Since  the  probing  nanorods  are  measured  in  nanometers,  this 
method  can  be  used  for  the  in  vivo  mapping  of  the  rheological  properties  of  biofluids  and  polymers  on  a  microscopic  level  at 
short  time  intervals  when  other  methods  fall  short. 


1.  INTRODUCTION 

Many  polymers  and  biofluids  rapidly  react  on  the  environ¬ 
mental  conditions  by  changing  their  rheological  properties.1 
Rheology  of  polymer  solutions/gels  (e.g.,  hazardous  fluids 
containing  thickeners)  and  biofluids  depends  on  the  concen¬ 
tration,  level  of  cross-linking,  oxygen  content,  temperature,  pH, 
and  many  other  environmental  parameters.  In  biofluids  such  as 
cellular,  mucosal,  and  tissue  fluids,  the  fluid  viscosity  changes 
over  time  making  it  challenging  to  quantify  and  analyze  these 
changes.  Drastic  thickening  of  fluids  over  a  short  time  interval 
often  results  in  far-reaching  consequences;  for  example,  sickle 
cell  disease  results  in  the  increase  of  cytosolic  viscosity  inside 
the  red  blood  cells,2,3  and  polymerization  of  fibrinogen  leads  to 
fibrin  clot  formation  and  wound  healing.4-6  In  polymer/ 
biomedical  engineering,  fast  polymerization  processes  are 
frequently  used  for  different  applications;  for  example,  by 
HEMA  (2-hydroxyethyl-methacrylate)  photopolymerization 
one  can  make  contact  disposable  lenses.7-  Alginate  and 
other  polysaccharides,  gelling  in  fractions  of  a  second,  are  used 
for  manufacturing  fibers  and  films.  These  broad  applications 
and  abundance  of  rapidly  gelling  systems  call  for  methods  of 
their  analyses  and  characterization.  The  physics  of  nucleation 
and  the  associated  rheological  fingerprints  of  the  gelation 
process  are  hidden  at  the  nanoscale  and  almost  indistinguish¬ 
able  at  the  microscale.  The  available  rheological  techniques  fall 
short  of  providing  viable  information  about  the  in  vivo  physico¬ 
chemical  reactions  and  gelation  mechanisms  in  the  real  time1 
especially  when  the  amount  of  liquids  is  limited  to  micro¬ 
droplets.3'12 
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A  number  of  the  microrheological  methods  were  recently 
developed  to  study  the  rheology  inside  the  living  cells.13-15 
These  methods  are  summarized  in  ref  16.  Typically,  a  probe  in 
the  form  of  a  solid  particle  is  embedded  into  the  cell  and  one 
measures  its  response  to  the  thermal  fluctuations  (passive 
microrheology15)  or  to  the  external  field,  for  example,  magnetic 
field  (active  microrheology14).  Active  microrheology  was  used 
to  study  the  rheological  properties  of  different  fluids.17  19 
Magnetic  rotational  spectroscopy  (MRS)  is  a  part  of  active 
microrheology  where  magnetic  nanorods  with  unique  rotational 
features  deserve  a  special  attention.20  26  In  the  past  decade,  the 
most  attention  has  been  paid  to  the  development  of  MRS  for 
simple  Newtonian  fluids.  In  this  Article,  we  develop  an  MRS 
which  can  be  used  to  study  the  rheological  properties  of  non- 
Newtonian  fluids  with  an  exponentially  fast-growing  viscosity. 
We  believe  that  the  MRS  described  here  will  provide 
biophysicists  and  chemists  with  a  unique  tool  to  probe  and 
map  the  physicochemical  properties  of  sophisticated  micro¬ 
systems  in  short  time  intervals. 

As  an  illustration  of  the  robustness  of  the  proposed  method, 
we  use  UV  polymerization  of  HEMA  to  analyze  its  time- 
dependent  viscosity.  The  HEMA  viscosity  exponentially 
increases  during  free  radical  (three-dimensional  UV)  polymer¬ 
ization  in  the  presence  of  the  cross-linker  DEGDMA 
(diethylene  glycol  dimethacylate).7  The  gel  point  for  the 
system  can  be  reached  in  60—200  s  depending  on  the 
DEGDMA  concentration.7  HEMA  polymerizes  linearly 
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through  the  carbon— carbon  double  bond  and  cross-links 
through  the  two  carbon— carbon  double  bonds  present  in 
DEGDMA  (see  Supporting  Information  and  Figure  SI  for 
details). 

The  time-dependent  viscosity  of  many  monomers  under¬ 
going  polymerization  is  typically  described  by  the  following 
equation:  ' 

>l(t )  =  %  et/r  (1) 

where  rj0  is  the  initial  viscosity  of  monomer  solution,  t  is  the 
time,  and  t  is  the  characteristic  time  of  polymerization.  This 
characteristic  time  is  considered  a  phenomenological  parameter 
of  the  polymerization  process.7’'1  The  physics  of  exponential 
thickening  of  fluid  via  three-dimensional  polymerization  is  not 
completely  understood,  yet  some  key  mechanisms  leading  to 
this  effect  have  been  identified.10 

In  order  to  assess  the  polymerization  kinetics  using 
rheological  data,  it  is  therefore  imperative  to  correlate  the 
kinetic  parameters  with  the  rheological  parameters  ?/0  and  t. 
The  existing  basic  MRS  model  does  not  take  into  account  the 
time  dependence  of  the  fluid  viscosity.  We  therefore  develop  a 
model  to  describe  the  dynamics  of  rotation  of  a  magnetic 
nanorod  in  a  fluid  with  an  exponentially  increasing  viscosity. 
Then,  we  apply  this  model  to  study  the  gelation  kinetics  of 
HEMA  during  photopolymerization.  We  analyze  the  exper¬ 
imental  data  on  rotation  of  magnetic  nanorods  and  perform  the 
FT-IR  measurements  of  the  HEMA  solutions  during  photo¬ 
polymerization.  All  these  experiments  show  the  robustness  of 
the  proposed  MRS  method  allowing  us  to  find  a  relation 
between  the  phenomenological  parameters  i/0  and  r,  concen¬ 
tration  of  the  cross-linker,  and  the  fraction  of  the  cross-linked 
material  at  different  stages  of  photopolymerization. 

2.  THEORY 

2.1.  Model.  The  movement  of  a  magnetic  nanorod  in  a 
rotating  magnetic  field  is  shown  in  Figure  1  (see  movie  S6  in 


Figure  1.  Counter  clockwise  rotation  of  nickel  nanorod  ( d  =  200  nm 
in  diameter)  in  a  plane.  Uniform  magnetic  field  rotates  with  constant 
frequency  f  =  0.5  Hz  forcing  the  nanorod  to  follow  it. 


Supporting  Information).  The  angle  a  specifies  the  direction  of 
applied  magnetic  field  B  with  respect  to  the  reference  axis  X. 
Since  the  applied  field  rotates  at  a  constant  rate,  this  angle 
depends  linearly  on  time,  a  =  2 nft,  where  f  is  the  frequency  of 
the  rotating  magnetic  field.  We  assume  that  the  magnetization 
vector  M  is  directed  parallel  to  the  nanorod  axis.  The  drag  force 
resists  the  nanorod  rotation  causing  the  nanorod  to  lag  behind 
the  field,  making  angle  0  with  vector  B  (Figure  l).  In  order  to 
derive  an  equation  governing  the  nanorod  rotation,  it  is 
convenient  to  count  its  revolutions  with  respect  to  the  fixed 
system  of  coordinates,  i.e.,  with  respect  to  the  X-axis.  Therefore, 
if  the  nanorod  axis  makes  the  angle  (p{t)  with  the  X-axis,  this 
angle  can  be  connected  with  a{t)  and  6{t)  as  q>{t)  =  In  ft  — 
0{t). 

The  torque  balance  equation  reads  (yd<}9/df)e  =  M  X  B, 

where  y  is  the  drag  coefficient,  and  e  is  the  unit  vector  directed 

perpendicularly  to  the  plane  of  the  nanorod  rotation.27,28 

Substituting  the  definition  of  angle  i p(t )  through  the  angles  a{t) 

and  0(t),  the  governing  equation  takes  on  the  following 

c  22.24 

form: 


mB  sin  0 


(2) 


where  m  is  the  magnetic  moment  of  the  nanorod.  The  drag 
coefficient  depends  on  the  nanorod  length  I,  its  diameter  d,  and 
liquid  viscosity  ?/(f)  =  t]0  exp (f/f)  as 


r(t)  = 


3  ln(//d)  -  A 


(3) 


where  constant  A  m  2.4  was  calculated  in  ref  29. 

In  the  case  of  a  constant  viscosity,  t  — *•  oo,  eq  2  has  a 
particular  solution  2 nfy  =  mB  sin  0.  This  solution  implies  that 
the  nanorod  follows  the  rotating  magnetic  field  B  making  angle 
0  with  the  field.  When  the  frequency  increases  above  the  critical 
frequency 

fa  =  mB  sin  0/ 2ny 

the  nanorod  slows  down  its  rotation.20,22  This  critical  frequency 
is  easy  to  observe  experimentally.  ’  '  Obtaining  this 
parameter  from  experiments,  one  can  infer  the  viscosity  of 
the  surrounding  fluid  by  solving  eq  4  for  y  and  then  applying  eq 
3  as  r  — >  oo  and  solving  it  for  >]0. 

When  the  viscosity  is  not  constant,  we  cannot  use  eq  4 
anymore.  Hence,  the  analysis  of  eq  2  deserves  special  care.  We 
analyze  this  case  by  introducing  the  dimensionless  times  T  = 
2 nft  and  T0  =  2nfr.  Then,  eq  2  can  be  rewritten  as  f)  eT7r°  ( 1  — 
A0/A T)  =  sin(0),  where 

n  _ 

mB[ 3  ln(//d)  —  A]  (5) 


After  making  this  nondimensionalization,  all  physical 
parameters  of  our  problem  collapse  into  two  dimensionless 
complexes,  f)  and  T0.  The  dimensionless  parameter  f)  describes 
all  possible  scenarios  of  the  nanorod  rotation:  assume  first  that 
the  drag  force  is  much  greater  than  the  magnetic  torque.  This 
implies  that  the  dimensionless  parameter  f)  is  much  greater 
than  1,  f  3=-  1.  Hence,  as  follows  from  the  governing  equation, 
(l  —  A0/AT)  =  (l /ff)e~T/r°  sin((9)  «  0,  or  0  »  T,  i.e.,  one 
would  not  expect  any  rotation  of  the  magnetic  nanorod.  In  the 
other  limit,  when  the  drag  force  is  much  smaller  than  the 
magnetic  torque,  the  /^-parameter  is  much  smaller  than  1,  /?  <3C 
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1.  As  follows  from  the  governing  equation,  /?  er/r°(l  —  d0/d T) 
=  sin(0),  one  would  expect  to  observe  the  nanorod  rotation 
until  the  left-hand  side  is  smaller  than  1. 

For  complete  analysis  of  the  nanorod  rotation,  it  is 
convenient  to  introduce  new  auxiliary  function  U  =  fSeT/T°. 
Then,  eq  2  is  rewritten  as  a  system  of  two  first-order  differential 
equations 

AO  I  AT  =  1  -  sin  0/U 

AU/AT  =  U/T0  (6) 


Dividing  the  first  equation  into  the  second,  we  obtain 


AB_ 

dU 


IT 


( 17  —  sin  0) 


(7) 


Equation  7  allows  one  to  perform  a  complete  analysis  of  the 
nanorod  dynamics  depending  on  the  viscosity,  nanorod 
magnetization,  magnetic  field,  and  initial  conditions.  According 
to  the  definition  of  function  17,  the  initial  condition  for  this 
function  reads  17( 0)  =  ft.  Therefore,  in  the  plane  (U,  0),  each 
vertical  line  U( 0)  =  /?  corresponds  to  a  continuous  set  of  initial 
angles  0(0)  =  0O.  Equation  7  has  many  solutions,  but  the 
nanorod  dynamics  is  described  by  a  single  solution 
corresponding  to  finite  interval. 

2.2.  Numerical  Analysis.  Dynamics  ofNanorods  Starting 
Rotation  at  cp(0)  =  0.  For  an  interpretation  of  the  images  of 
rotating  nanorods,  we  need  to  follow  the  <p(f)-dependence. 
Therefore,  we  rewrite  the  torque  balance  equation  in  the 
following  form  y(f)  Arp /At  =  mB  sin(2 jift  —  rp).  As  an  initial 
condition,  we  can  always  choose  cp( 0)  =  0  implying  that  the 
nanorod  is  initially  aligned  with  the  field.  The  right-hand  side  of 
this  equation  is  positive  within  the  time  interval  t  <  ta  where  tc 
is  the  first  root  of  the  equation  2 njtc  —  cp{tc )  =  0.  This  root  is 
well-defined  for  each  r:  it  corresponds  to  the  first  maximum  of 
function  rp{t ).  In  Figure  2a,  we  show  the  behavior  of  the 
nanorod  with  length  l  =  5  pm,  diameter  d  =  0.2  pm,  initial 
viscosity  of  the  HEMA  solution  ;/0  =  2.7  mPa-s,  magnetic  field  B 
=  0.0015  T,  and  frequency  of  the  rotating  field /=  0.5  Hz. 

In  Figure  2a,  the  nanorod  revolves  seven  times  at  the 
frequency  of  the  rotating  field,  and  the  dashed  and  solid  lines 
corresponding  to  the  field  and  nanorod  rotation  angles, 
respectively,  coincide.  On  the  other  hand,  the  dash— dot  line 
corresponding  to  the  0-angle  remains  at  zero.  Thus,  during  this 
rotation  period,  the  field  and  magnetization  are  oriented 
parallel  to  each  other.  After  seven  turns,  the  slope  of  the  red 
line  changes,  implying  a  gradual  decrease  of  the  rotation  rate 
from  0.5  to  0.45  Hz.  It  can  be  seen  from  the  graphs  that  the 
nanorod  ( rp )  follows  the  rotation  of  the  field  (a)  for  about  JVcr 
=  11.17  x  11  revolutions.  After  Ncr  turns,  the  derivative  of 
( p(N )  changes  sign  from  positive  to  negative  and  the  nanorod 
sways  back  for  about  45°  (see  movie  S7  in  Supporting 
Information).  After  that,  the  nanorod  rotates  back  and  forth 
with  decreasing  amplitude  and  then  completely  stops  at  N  x 
15.  In  experiments,  NCI  can  be  found  by  measuring  the  time  (or 
number  of  turns)  that  the  nanorod  takes  on  before  it  sways 
back  for  the  first  time.  In  Figure  2b,  we  plot  the  dependence  of 
this  critical  parameter  Ncr  as  a  function  of  T0  and  f). 

One  can  infer  from  these  dependencies  that  the  number  of 
turns  significantly  decreases  as  the  initial  viscosity  increases  and 
as  the  external  field  increases.  This  analysis  suggests  an 
experimental  strategy  for  the  investigation  of  rheological 
properties  of  polymeric  fluids  during  polymerization.  That  is, 


P 


Figure  2.  (a)  Numerical  solutions  for  rp,  a,  and  9  as  functions  of  the 
number  of  turns,  N  =  tf,  of  a  nanorod.  The  inset  shows  oscillating 
behavior  of  a  nanorod  which  is  about  to  stop  its  rotation  after  Ncr  ~  1 1 
revolutions,  (b)  Dependence  of  the  critical  number  of  turns  of  a 
nanorod  on  parameter  /). 


one  needs  to  measure  Ncr  and  use  t  as  an  adjustable  parameter 
to  fit  the  experimental  Ncr  with  the  theoretical  one. 

Dynamics  of  Nanorods  Starting  Rotation  at  Arbitrary 
Angle  0O.  There  are  three  possible  scenarios  of  nanorod 
rotation  starting  at  an  arbitrary  angle  90.  The  results  of 
numerical  analysis  of  eq  7  are  summarized  in  the  form  of  the 
phase  portrait  shown  in  Figure  3.  The  initial  conditions 
satisfying  the  inequality  U0  =  f)  <  1  correspond  to  the  case 
when  the  magnetic  torque  is  stronger  than  that  caused  by  the 
viscous  drag.  The  shaded  region  under  curve  U  =  sin  0 
corresponds  to  the  conditions  that  cause  the  rod  to  sway 
toward  the  field  direction  at  the  initial  instance  of  time.  When 
the  viscous  drag  takes  over,  the  rod  slows  down  its  rotation  and 
eventually  stops.  Since  magnetic  field  keeps  revolving,  the  angle 
0(17)  increases  with  each  revolution  even  if  the  nanorod  is  not 
moving.  Therefore,  the  integral  curves  coming  out  from  the 
shaded  region  first  decrease,  form  minimums  at  curve  U  =  sin  0, 
and  then  turn  up  to  increase. 

This  behavior  can  be  easily  observed  by  first  aligning  the 
nanorods  with  a  bias  field  and  then  applying  the  rotating  field 
which  is  directed  perpendicular  to  the  nanorods  at  t  =  0. 
Therefore,  the  initial  condition  for  this  experiment  is  stated  as 
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Figure  3.  Phase  portrait  for  eq  6  showing  the  solution  behavior  for 
different  initial  conditions  170  and  00  at  fixed  dimensionless  T0  =  n. 
Only  17  >  0  is  of  practical  interest;  hence,  the  integral  curves  for  U  <  0 
are  not  shown. 


Oq  =  n/'l.  The  time  when  the  nanorods  form  the  minimum 
angle  #min  with  the  rotating  field  can  be  thus  obtained  as 
jjgTmmjTo  =  sin  g/m.n  or  rmin  =  2nfc  ln(sin  0mJ/3).  Since  0mm  and 
Tmin  are  well-defined  experimentally,  one  can  infer  the  z 
parameter  from  this  relation. 

Numerical  analysis  shows  that  all  curves  coming  out  of  the 
shaded  region  first  pass  a  minimum,  then  find  their  way  to 
move  at  a  close  proximity  of  curve  (A).  When  17  is  large,  U 
1,  eq  7  can  be  integrated  to  give  0  =  T0  In  17,  which  means  that 
angle  0  increases  linearly  with  time.  These  asymptotes  explains 
why  all  integral  curves  run  parallel  to  each  other  as  U  increases. 

The  integral  curves  coming  out  from  the  regions  outside  the 
shaded  area  in  Figure  3,  always  increase  as  17  increases. 
Therefore,  as  follows  from  the  definition  of  angle  0,  we  have  the 
following  inequality  d 0/ d t  =  2 nj  —  dcp/ dt  >  0.  Rewriting  this 
inequality  for  the  observable  angle,  we  obtain  2 nj>  d(p/dt.  This 
behavior  suggests  that  the  external  field  rotates  faster  than  the 
rod  for  these  initial  conditions.  The  phase  portrait  sketched  in 
Figure  3  is  applicable  for  all  T0,  yet  the  slopes  of  the  integral 
curves  are  influenced  by  this  parameter  T0. 

The  characteristic  features  of  the  nanorod  behavior  can  be 
further  elucidated  by  following  the  trajectories  belonging  to 


different  regions  on  the  phase  portrait.  For  example,  the  initial 
conditions  0O  =  0,  2,  and  2.9  correspond  to  three  different  types 
of  the  nanorod  trajectories.  From  equation  (p(t )  =  a(t )  —  0(t), 
one  can  infer  the  following  initial  condition  for  ^-solutions: 
(p(0)  =  —0(0).  The  results  of  these  calculations  are  presented  in 
Figure  4  for  physical  parameter  z  =  4. 

Figure  4a  shows  three  different  dependencies  of  angle  0  on 
parameter  17  at  different  initial  conditions  0O,  and  the  same 
L7(0)  =  0.5.  In  region  U  >  sin(0o),  the  initial  viscous  drag  is 
stronger  than  the  magnetic  torque;  hence,  the  field  revolves 
faster  than  the  nanorod.  As  follows  from  eq  7,  the  derivative 
d0/dU  is  positive.  With  time,  the  viscous  drag  increases 
exponentially;  hence,  the  field  keeps  revolving  faster  than  the 
nanorod  rotates.  For  initial  condition  0O  =  0,  the  magnetic  field 
and  nanorod  magnetization  are  initially  coaligned  (Figure 
4a,b,c).  As  shown  in  Figure  4b,  after  approximately  0.2 
revolutions,  the  angle  0  stops  changing.  This  implies  that  the 
nanorods  are  caught  by  the  field  for  some  time  (Figure  4c). 
Then,  after  rotation  in  unison  with  the  field  for  approximately  2 
revolutions,  the  frequency  of  nanorod  rotation  starts  to  slow 
down.  The  nanorod  sways  back  for  the  first  time  at  NCI  =  4. 

The  dash— dot  trajectory  corresponding  to  the  initial 
condition  0O  =  2  in  Figure  4  a,b,c  shows  the  distinct  features 
of  the  trajectories  with  initial  conditions  taken  from  region  U  < 
sin(#0).  In  this  region,  the  magnetic  torque  is  stronger  than  the 
viscous  drag.  Rotating  magnetic  field  attempts  to  turn  the 
nanorod  toward  the  field  direction.  Simultaneously,  the  field 
keeps  rotating  toward  the  nanorod.  This  results  in  a  decrease  of 
the  angle  0  as  dictated  by  eq  7  showing  negative  derivative  d0/ 
dU  in  the  region  in  question.  After  approximately  0.2 
revolutions,  the  angle  0  stops  changing  (Figure  4b).  This 
means  that  the  nanorods  are  caught  by  the  field  for  some  time 
repeating  the  behavior  of  nanorods  which  started  at  0O  =  0  (the 
dash— dot  line  in  Figure  4c).  Then,  after  rotation  in  unison  with 
the  field  for  approximately  2  revolutions,  the  frequency  of 
nanorod  rotation  starts  to  slow  down.  The  nanorod  sways  back 
for  the  first  time  at  Na  =  4. 

The  dotted  trajectory  corresponding  to  the  initial  condition 
0o  =  2-9  in  Figure  4  a,b,c  elucidates  the  distinct  features  of  the 
trajectories  with  the  initial  conditions  taken  from  region  17  > 
sin(#0).  In  this  region,  the  drag  force  is  greater  than  the 
magnetic  torque.  The  explanation  of  the  nanorod  behavior  is 
analogous  to  that  discussed  above  for  17  >  sin(0o),  0O  =  0. 


Figure  4.  Nanorod  behavior  at  three  initial  conditions,  60  =  0,  2,  and  2.9  radians,  (a)  Phase  portrait  of  eq  7  showing  the  solution  behavior  for  these 
initial  conditions,  (h)  Numerical  0-solutions  as  a  function  of  the  number  of  turns,  N  =  t-f.  (c)  Numerical  (^-solutions  as  a  function  of  the  number  of 
turns. 
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Figure  5.  Nanorod  behavior  at  different  initial  conditions  /?  =  U( 0)  =  0.3,  0.5,  0.7.  (a)  Phase  portrait  of  eq  7  showing  the  solution  behavior  for 
different  initial  conditions  17( 0)  at  60  =  0.  (b)  Numerical  ^-solutions  as  a  function  of  the  number  of  turns,  N  =  t-f.  (c)  Numerical  (^-solutions  as  a 
function  of  the  number  of  turns. 


For  all  these  three  cases,  the  parameter  Ncr  is  the  same,  and  it 
is  equal  to  Ncr  =  4.  From  the  numerical  analysis,  we  conclude 
that  the  behavior  of  nanorods  at  the  first  moments  of  time  does 
depend  on  angle  60,  but  the  critical  number  of  turns  NCI  does 
not  depend  on  the  initial  orientation  of  the  nanorods. 

Effect  of  [/-Parameter  on  Nanorod  Dynamics.  Effect  of  the 
/^-parameter,  i.e.,  the  initial  condition  U( 0),  can  be  explained  by 
following  the  trajectories  corresponding  to  ft  =  U0  =  0.3,  0.5, 
and  0.7  and  assuming  that  the  initial  angle  60  =  0.  The 
dependencies  of  functions  6  and  q)  on  time  are  presented  in 
Figure  5  where  we  have  chosen  t  =  4. 

Figure  5  reveals  that  the  nanorod  behavior  for  different  (3 
does  not  change  significantly:  the  behavior  of  the  ^-solutions 
on  the  number  of  turns  is  very  similar  among  the  solutions.  The 
same  is  valid  for  the  ^-solutions.  At  initial  moments  of  time,  the 
q)-  and  cr-solutions  coincide,  i.e.,  the  nanorods  are  coaligned 
with  the  field  and  revolve  with  the  same  frequency  as  that  of  the 
driving  magnetic  field  (Figure  5c).  The  nanorod  keeps  rotating 
until  it  reaches  a  critical  number  of  turns,  Ncr  .  This  critical 
number  depends  on  the  initial  conditions;  for  example,  Ncr( 0.3) 
=  2.8,  Ncr(0.5)  =  4,  Ncr( 0.7)  =  6.1.  The  lower  the  /^-parameter, 
the  greater  the  critical  number  of  turns.  As  follows  from  eq  5, 
the  /^-parameter  can  be  decreased  by  decreasing  initial  viscosity, 
or  by  making  the  nanorod  shorter,  or  by  applying  a  stronger 
magnetic  field,  or  by  decreasing  the  frequency  of  magnetic  field 
rotation. 

This  numerical  analysis  of  eq  7  allows  one  to  design  the  MRS 
experiments  with  different  materials  and  estimate  the  critical 
parameters  such  as  the  critical  number  of  turns  or  r-parameter 
specific  for  different  liquid— nanorod  pairs. 

3.  MATERIALS  AND  METHODS 

HEMA  monomer  (Sigma- Aldrich  Co.  LLC)  was  mixed  with  the  cross- 
linking  agent  DEGDMA  (Sigma- Aldrich  Co.  LLC).  Five  solutions  with 
different  concentrations  of  the  cross-linking  agent,  1.86,  3,  4.5,  7,  and 
9.3  wt  %,  were  prepared  for  testing.  The  photoinitiator  2,2-dimethoxy- 
2-phenylacetophenone  (Sigma- Aldrich  Co.  LLC)  of  1  wt  %  relative  to 
the  monomers  was  added  to  the  solutions.  Nickel  nanorods  were 
synthesized  inside  pores  of  alumina  membranes  (Watman  Ltd.)  by  the 
electrodeposition  technique.23'31  The  synthesis  is  described  in  detail  in 
Supporting  Information.  HEMA  solution  was  added  to  the  beaker  with 
nickel  nanorods,  and  an  ultrasound  agitation  was  used  to  prepare  a 
uniform  suspension  of  nanorods.  Concentration  of  nanorods  in  the 
suspension  was  0.05  wt  %.  In  order  to  exclude  the  interactions 
between  the  nanorods,  we  always  double-checked  that  the  chosen 
nanorod  was  far  away  from  its  neighboors.  The  distance  to  the  closest 
nanorod  was  at  least  7  times  greater  than  the  nanorod  length  (see 


Figure  S9  in  Supporting  Information).  This  guarantees  a  weak 
dependence  of  nanorod  behavior  on  dipole— dipole  interactions  as 
discussed  in  ref  32.  A  1.8  pL  droplet  of  the  nanorod  dispersion  was 
placed  between  two  glass  slides.  A  64-^m-thick  Kapton  tape  was 
inserted  between  the  ends  of  two  glass  slides  to  keep  them  separated. 
The  drop  was  sandwiched  between  glass  slides  and  was  investigated 
under  the  50X  objective  of  the  BX-51  Olympus  microscope.  We 
applied  the  dark  field  microscopy  technique  to  detect  the  200-nm- 
diameter  nanorods  which  are  not  visible  in  the  bright  field.33  First,  the 
microscope  was  focused  on  the  upper  glass  slide,  and  then  by  moving 
the  objective  down,  it  was  focused  on  the  bottom  glass  slide.  This  way, 
it  was  possible  to  focus  the  objective  on  the  magnetic  nanorod  which 
was  situated  in  the  middle  of  the  optical  cell.  Each  experiment  was 
started  by  confirming  that  the  tested  nanorod  was  situated  at  least  20 
/ im  away  from  each  optical  slide  and  the  wall  effects  can  be  safely 
neglected.  This  nanorod  position  did  not  change  during  the  HEMA 
photopolymerization  experiment,  i.e.,  the  nanorods  have  never  settled 
down.  The  following  estimates  explain  these  observations.  Assuming 
that  the  nanorod  is  settling  down  at  velocity  v  and  it  keeps  the  long 
axis  parallel  to  the  plates,  the  drag  force  on  a  nanorod  is  written  as  F  = 
47njlv/\n(l/d).34  This  drag  force  is  balanced  by  the  buoyancy  force,  F  = 
ndPlApg/ 4,  where  A p  =  8000  kg/ m3  is  the  density  difference  and  g  is 
acceleration  due  to  gravity.  Therefore,  the  settling  velocity  is  estimated 
as  v  =  d2lApg  \n(l/d)/(l6r]l)  =  d2Apg  \n(l/d)/(l6rj).  Assuming  l/d  = 
30,  d  =  0.2  pm,  and  fj  =  2  mPa-s  (viscosity  prior  to  photo 
polymerization),  we  have  v  =  0.33  pm/s.  If  at  the  first  moment  a 
nanorod  was  positioned  in  the  middle  of  the  optical  cell,  which  is  30 
pm  away  from  the  bottom  plate,  it  would  take  about  t  =  30/0.33  =  91 
s  to  settle  down.  As  shown  in  the  experiments,  the  HEMA 
photopolymerization  was  initiated  in  less  than  10  s  after  placing  the 
droplet  on  the  glass  slide,  and  it  took  about  40  s  to  change  the  HEMA 
viscosity  100  times.  This  explains  why  the  buoyancy  effect  was 
insignificant  and  why  we  were  able  to  follow  the  same  nanorod  during 
photopolymerization.  A  supplementary  video  (Supporting  Information 
S8)  shows  the  behavior  of  a  nanorod  subject  to  the  Brownian 
fluctuations  during  HEMA  photopolymerization  in  the  absence  of 
magnetic  field.  Mostly  in-plane  fluctuations  were  observed,  and  the 
nanorod  persistently  stayed  in  focus  while  the  film  was  recorded. 

Two  magnetic  coils  orthogonal  to  each  other  were  used  to  produce 
a  rotating  magnetic  field.  We  used  E-66-100  coils  (15  mm  diameter 
and  8.5  mm  thick,  Magnetic  Sensor  Systems).  The  current  was 
controlled  with  a  LabView  program.  The  design  of  the  setup  is 
described  in  ref  23.  The  magnetic  field  at  the  center  of  the  optical  cell 
was  estimated  to  be  B  =  1.5  X  10-3  T.  A  SPOT  videocamera  (SPOT 
Imaging  Solutions,  Inc.)  was  used  for  filming  the  rotating  nanorod 
during  HEMA  polymerizaion.  In  order  to  initiate  the  polymerization, 
the  B-100SP  UV  lamp  (UVP,  LLC)  was  placed  100  mm  apart  from  the 
glass  slides  containing  the  HEMA  solution.  The  light  intensity  at  this 
distance  was  approximately  10  mW/cm2.  The  captured  videos  were 
analyzed  using  VirtualDub  software  (http://www.virtualdub.org). 


E 


dx.doi.org/l  0.1 021  /Ia301 9474 1  Langmuir  XXXX,  XXX,  XXX-XXX 


Article 


Langmuir 


Magnetic  properties  of  nickel  nanorods  vary  significantly  depending 
on  the  parameters  of  electrodeposition  process.35'36  In  order  to 
calculate  the  magnetic  moments  of  the  nickel  nanorods,  we  conducted 
the  MRS  experiments  in  ethylene  glycol,  a  simple  Newtonian  fluid 
with  )/  =  16  mPa'S,  and  used  eq  4  to  calculate  m  (for  details,  see 
Supporting  Information).  The  obtained  magnetic  moment  appeared  to 
be  approximated  well  by  the  formula  m  =  VB///J0  where  V  is  the 
nanorod  volume,  /v0  is  the  magnetic  permeability  of  vacuum,  and  /  is  a 
constant,  /  =  92.  We  also  calibrated  the  MRS  viscosity  measurements 
on  fluids  with  known  viscosities.  First,  we  measured  the  bulk  viscosity 
of  glycerin  ( I /  =  778  mPa-s)  and  its  water  mixtures  using  DV-III  Ultra 
Rheometer  (Brookfield  Eng.)  and  then  compared  the  data  with  the 
MRS  measurements.  The  results  are  in  a  good  agreement  (see  S3  in 
Supporting  Information).  These  experiments  confirmed  that  we  can 
measure  viscosities  up  to  780  mPa-s  with  our  setup. 

After  calibration  of  our  MRS  system,  we  obtained  the  characteristic 
viscosity  i]0  by  rotating  the  nanorod  in  HEMA  solution  prior  to 
exposing  it  to  the  UV  radiation.  Once  the  critical  frequency  was 
reached,  we  solved  eq  4  for  the  viscosity  to  obtain  >]0  =  2.7  mPa-s.  This 
characteristic  viscosity  did  not  change  when  we  added  the  cross-linker. 

In  all  MRS  experiments,  we  were  trying  to  measure  the  phase  lag 
between  the  applied  signal  and  nanorod  displacement.  Our  system  was 
unable  to  detect  any  phase  lag  until  the  nanorod  reached  NCI 
revolutions.  Therefore,  within  this  time  interval  the  material  reacted 
as  a  purely  viscous  liquid,  and  no  fingerprints  of  material  elasticity  were 
observed.  With  this  confirmation  in  hand,  we  applied  our  model  to 
interpret  the  experimental  data. 

4.  RESULTS 

With  the  obtained  parameters,  we  used  the  following  strategy 
to  evaluate  the  dependence  of  characteristic  time  r  on  the 
cross-linker  concentration  during  UV-polymerization.  For  each 
of  five  solutions,  corresponding  to  1.86,  3,  4.5,  7,  and  9.3  wt  % 
cross-linker  concentrations,  we  obtained  the  critical  number  of 
nanorod  revolutions  Ncrexp.  Ten  different  nanorods  with  the 
lengths  ranging  from  3.9  to  11.85  /(m  were  used  in  these 
experiments  for  each  droplet.  The  diameter  of  all  nanorods  was 
controlled  by  the  diameter  of  pores  in  alumina  membrane  and 
it  was  equal  to  200.57  ±  2.59  nm  (Figure  S2(c)  in  Supporting 
Information).  Then,  we  ran  the  numerical  analysis  of  our 
model,  eqs  6,  to  fit  the  experimental  data  using  t  as  an 
adjustable  parameter  of  the  model.  The  obtained  T-values  were 
used  to  plot  ten  exponential  curves  for  each  cross-linker 
concentration. 

In  Figure  6a,  we  plot  relative  viscosity  f](t)/rj0  of  HEMA  with 
1.86  wt  %  concentration  of  the  cross-linker  measured  with  the 
MRS  and  compared  these  measurements  with  available  data 
reported  in  ref  7  for  the  different  powers  of  UV  irradiation.  It 
can  be  seen  that  the  MRS  results  are  well-fitted  by  the 
exponential  curve  and  the  relative  viscosity  curve  for  our  UV 
source  (W  =  10  mW/cm2)  is  situated  between  those  obtained 
in  ref  7  at  stronger  UV  irradiation  with  W=  40  mW/cm"  power 
and  weaker  UV  irradiation  with  W  =  0.25  mW/cm2  power. 
Therefore,  the  obtained  results  appear  reasonable  and  coherent 
with  the  available  data. 

In  Figure  6b,  we  plot  five  dependencies  of  the  relative 
viscosity,  each  being  related  to  the  given  cross-linker  content. 
Any  point  on  these  curves  represents  the  average  value  taken 
over  ten  measurements  for  each  concentration  of  the  cross¬ 
linker.  Error  bars  provide  the  standard  deviation  from  the  mean 
value.  In  the  inset,  we  show  the  r-dependence  on  the  cross¬ 
linker  concentration. 

It  can  be  seen  from  the  graphs  that  the  fastest  change  of 
viscosity  is  observed  in  the  system  with  the  highest 
concentration  of  the  cross-linker.  The  dependence  of  the 
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Figure  6.  (a)  Viscosity  change  during  UV  polymerization  of  HEMA 
with  1.86  wt  %  concentration  of  the  cross-linker  measured  with  MRS 
(/  =  0.5  Hz,  r  =  17  s)  compared  with  bulk  measurements  of  ref  7 
corresponding  to  r  =  5.8  s  and  r  =  34.5  s.  (b)  Viscosity  change  during 
UV  polymerization  for  different  concentrations  of  the  cross-linker 
(from  the  right  to  the  left  1.86,  3,  4.5,  7,  9.3  wt  %).  Inset:  the 
dependence  of  characteristic  time  t  on  the  weight  fraction  of  the  cross¬ 
linker,  C,  in  the  monomer/cross-linker  mixture  (diamonds  are  the 
experimental  points);  line  is  a  cubic  polynomial  fit.  (c)  Theoretical 
dependence  of  the  critical  number  of  turns  Ncr  of  a  nanorod  as  a 
function  of  characteristic  time  t. 
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characteristic  time  r  on  the  concentration  of  the  cross-linker 
can  be  fit  with  a  cubic  polynomial  r  =  a-C-1/3,  where  C  =  mcr/ 
mmix  is  the  weight  fraction  of  the  cross-linker  of  weight  mcr  in 
the  monomer/cross-linker  mixture  of  weight  mmix  and  a  =  4.36 
s  is  a  coefficient  determined  from  the  fit.  A  theoretical 
dependence  of  Na  on  r  is  shown  in  the  Figure  6c. 

This  analysis  demonstrates  that  the  magnetic  rotational 
spectroscopy  allows  one  to  measure  the  viscosity  of 
polymerized  solution  when  its  viscosity  drastically  changes 
over  a  short  time  interval.  However,  from  the  MRS  data,  one 
cannot  resolve  whether  the  polymer  has  already  undergone  the 
transition  to  a  gel  form  or  it  is  still  a  liquid  and  not  completely 
cross-linked  yet.  We  therefore  conducted  a  series  of  additional 
experiments  enabling  one  to  distinguish  the  conformational 
changes  of  HEMA  solutions  at  different  moments  in  time. 

5.  DISCUSSION 

HEMA  polymerizes  through  the  carbon— carbon  double  bond 
and  cross-links  through  the  two  double  bonds  in  DEGDMA 
(see  Figure  SI  in  Supporting  Information).  Therefore,  in  order 
to  evaluate  the  cross-linking  kinetics,  we  used  FTIR  spectros¬ 
copy  (Nicolet  Magna  550).  We  followed  the  rate  of  decrease  of 
carbon— carbon  double  bonds  in  the  system  by  measuring  the 
rate  of  disappearance  of  the  1635  cm-1  peak  corresponding  to 
the  carbon— carbon  double  bonds.37  The  intensities  of  the 
absorbance  peaks  were  measured  before  polymerization  (solid 
line  in  the  inset  of  Figure  7)  and  after  30,  50,  and  60  s  of  the 
UV  exposure.  The  dashed  line  in  the  inset  of  Figure  7  shows 
the  intensity  of  the  1635  cm-1  peak  after  60  s  of  photo¬ 
polymerization. 
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Figure  7.  Solid  line  illustrates  the  change  of  relative  viscosity  with  time 
of  photopolymerization  of  the  4.5  wt  %  cross-linker  solution  (left  y- 
axis),  the  stars  show  the  conversion  of  double  bonds  during 
photopolymerization  (right  v-axis),  and  the  dashed  line  is  a  trend 
line.  The  inset  shows  a  FTIR  spectrum  of  the  solution  near  1635  cm-1 
before  (solid  line)  and  after  60  s  of  polymerization  (dashed  line). 


The  1720  cm-1  peak  corresponding  to  carbonyl  group  in 
HEMA  and  DEGDMA  was  used  as  a  reference  to  calculate  the 
conversion  of  carbon— carbon  bonds  in  the  system,  because  it 
does  not  change  during  the  photopolymerization.  Details  on 
the  calculations  from  the  FTIR  measurements  can  be  found  in 
Supporting  Information  (see  Figure  S4).  We  followed  the 
double  bond  conversion,  (f>,  corresponding  to  the  ratio  of  the 
reacted  double  bonds  to  the  total  number  of  double  bonds  in 
the  system.  The  conversion  was  found  from  the  following 


equation37  (f>  =  100  X  (l  —  XtS0/X0St),  where  X0  and  Xt  are  the 
areas  under  the  1638  cm-1  peak  at  time  0  and  time  t, 
respectively(integration  interval  is  1604— 1656  cm-1),  and  S0 
and  St  are  the  areas  under  the  reference  1720  cm-1  peak  at  time 
0  and  time  t,  respectively  (integration  interval  is  1652—1774 
cm-1).  The  conversion  of  the  double  bonds  (p  with  time  is 
plotted  in  Figure  7  as  the  dashed  line.  We  found  that  the 
number  of  double  bonds  in  the  solution  has  decreased  by 
91.4%  after  60  s  of  photopolymerization.  This  reduction  of  the 
total  number  of  the  double  bonds  in  the  system  is  caused  by 
both  the  linear  polymerization  of  HEMA  and  cross-linking  of 
HEMA  through  the  double  bonds  of  DEGDMA.  However,  the 
FTIR  analysis  (showing  only  the  total  extent  of  the 
polymerization)  does  not  provide  information  on  the  amount 
of  the  polymeric  material  involved  in  the  cross-linked  gel. 
Therefore,  we  determined  the  amount  of  gel  fraction  in  the 
polymerizing  system  as  a  function  of  time  (Table  l).  This 


Table  1.  Gel  Fraction  at  Different  Stages  of  the 
Photopolymerization 

time  of  the  UV  exposure,  s  fraction  of  cross-linked  polymer  at  time  t,  wt  % 
45  (tcrl  in  Figure  7)  23.1  ±  2.11 

50  (tcr2  in  Figure  7)  30.9  ±  2.8 

55  (1ct3  in  Figure  7)  35.7  ±1.4 


fractionation  was  achieved  by  washing  out  the  cross-linker, 
photoinitiator,  and  a  linear  polymer  from  the  sample 
immediately  after  the  photopolymerization.  The  experimental 
protocol  is  described  in  Supporting  Information. 

Three  points  in  Figure  7,  taken  at  tcrl,  fcr2,  and  tcr3  correspond 
to  the  time  moments  when  nanorods  of  different  lengths  sway 
back  for  the  first  time.  As  follows  from  the  analysis  of  Figure  6c, 
the  longer  nanorods  having  greater  parameter  stop  rotating 
earlier  than  the  shorter  ones.  Accordingly,  the  6.5  Jim  long 
nanorods  stopped  at  time  fcrl,  the  5.8  /im  long  nanorods 
stopped  at  time  tcr2,  and  then  the  3.9  f-im  long  nanorods 
stopped  at  time  fcr3.  The  data  presented  in  Table  1  suggest  that 
the  longer  nanorods  most  likely  probe  HEMA  having  weakly 
cross-linked  gel  nuclei,  while  the  shorter  nanorods  are  able  to 
probe  HEMA  having  stronger  cross-linked  gel  domains  with 
the  total  amount  of  gel  in  the  suspension  reaching  ~36  wt  %. 
Thus,  the  MRS  provides  instructive  information  about  the 
viscosity  of  polymers  undergoing  gel  formation. 

6.  CONCLUSIONS 

We  theoretically  developed  and  implemented  a  new  MRS 
method  which  can  be  used  for  the  in  situ  (or  in  vivo) 
rheological  measurements  of  biofluids  and  polymer  systems 
when  the  fluid  viscosity  increases  with  time  exponentially  fast. 
Only  a  small  quantity  of  the  sample  is  needed:  for  example, 
assuming  that  the  rotating  5-^m-long  nanorod  can  cover  the  1 
pL  volume,  one  can  measure  the  rheological  property  of  a  drop 
of  comparable  size.  The  material  rheology  can  be  probed  on 
micro-  as  well  as  on  nanolevels,  depending  on  the  size  of  the 
nanorods  used.  Remarkably,  an  exponential  increase  of  viscosity 
can  be  traced  beyond  the  point  when  the  polymer  system 
undergoes  the  transition  to  a  gel  and  the  domains  start  to 
appear.  We  showed  that  the  MRS  with  nanorods  probes  the 
viscosity  of  HEMA  when  the  conversion  of  double  bonds 
reaches  ~90%.  We  expect  that  this  method  will  open  new 
horizons  in  the  quantitative  rheological  analysis  of  fluids  inside 
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the  living  cells,  microorganisms,  and  aerosol  droplets  with 
thickeners. 
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